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Abstract 




In this paper, we reexamine the new physics contributions to seven well 
measured B — > PP decays in the standard model (SM) and the general two- 
Higgs-doublet model (model III) and compare the theoretical predictions with 
the new data. Within the considered parameter space we found that: (a) the 
measurements of the branching ratios for B — ► -k + it~ , K~ir + and K°ir + lead 
to a strong constraint on the form factor Fq w (0): Fo n (0) = 0.24 ± 0.03; (b) 
the new physics enhancements to the penguin-dominated B — > Kir and Krj' 
decays are significant in size, ~ (40 — 65)% w.r.t the SM predictions, and 
play an important role in restoring the consistency between the data and the 
theory. 
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As is well known, one of the main objectives of B experiments is to probe for the possible 
effects of new physics beyond the standard model (SM). Precision measurements of the B 
meson system can provide an insight into very high energy scales via the indirect loop effects 
of the new physics |]J^] • 

Up to now, CLEO, BaBar and Belle Collaboration [Q-|5|] have observed eighteen two-body 
charmless hadronic B u ^ meson decay modes. Seven well measured B — > PP ( here P refers 
to the light pseudo-scalar mesons ) decays are 



B — > 7r ± 7r T , Kir, Krj 



(1) 



These decay modes are closely related through the isospin symmetry and the £77(3) flavor 
symmetry, phenomenologically very interesting due to their key role in extracting the unitary 
angles a and 7, and the appearance of the so-called rj'K puzzle: the observed B — ► Kr/' 
decay rates |3|-^| are much larger than what expected in the standard model based on the 
effective Hamiltonian with generalized factorization (GF) approach To accommodate 

the data, one may need an additional contribution unique to the rf meson in the framework 
of the SM [p7J| -[T3]| , or enhancements from new physics models beyond the SM [ I4|,|15 . 



In a previous paper JT3J], we considered the second possibility and calculated the new 
physics effects on the two-body charmless hadronic B meson decays in the general two-Higgs- 
doublet models (2HDM's) ||16|| . In this paper, we focus on seven well measured B — > PP 
decay modes and compare the theoretical predictions with the newest data. We still use 
the low-energy effective Hamiltonian with GF approach |8],|9| to calculate the new physics 
contributions. For recent studies of B — > PP decays in the SM with the QCD factorization 
(BBNS) and perturbative QCD (PQCD) approaches Hl7| , [18| , one can see the papers [|Kf -|22" 
and references therein. 

For the inclusive three-body decays b — > sqq with q e {u, d, s} the effective Hamiltonian 
can be written as lELHI, 



Gp I 



H eff (AB = 1) = ^ ^£0, [V ub V* s Q] + V cb V*Q 



10 



^ CjQj + CgQ & 



i=3 



• (2) 



The explicit expressions for all operators can be found easily for example in Ref. [§]. For 
b — > dqq decays, one simply makes the replacement s — > d. Following Ref. ||, we also neglect 
the effects of the electromagnetic penguin operator Qi^, the weak annihilation and exchange 
diagrams. Within the SM and at scale M w , the Wilson coefficients Ci(M w ), • • • , C W (M W ) 
at next-to-leading logarithmic order NLO) and C g (M\y) at leading logarithmic order (LO) 
have been given for example in Ref. . 

In a recent paper ||23j| , Chao et al. studied the decay b — » 57 in model III (the third 
type of 2HDM's) by assuming that only the couplings X u = \Xtt\e tdt and A{,& = lA^le^ 6 are 
non-zero. They found that the constraint on M#+ imposed by the CLEO data of b —>■ 57 
can be greatly relaxed by considering the phase effects of \ tt and A&&. From the studies of 
Refs. P,p3]|, we know that for model III the parameter space 



A 



v - 
A«| 



0, for ij 7^ tt, or bb, 
= 0.3, |A 66 | = 35, 6 = (0° 



30 c 



M H+ = (200 ± 100) GeV, 



(3) 
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are allowed by the available data, where 6 = 6bb — Ott- In this paper, we calculate the new 
physics contributions to seven B meson decay modes in the Chao-Cheung-Keung (CCK) 
scenario of model III [[23[ . 

Following the same procedure as in the SM, it is straightforward to calculate the new 
7-, Z°- and gluonic penguin diagrams induced by the exchanges of charged-Higgs bosons 



appeared in model III (for details of the calculations, see Ref. [HI). In the NDR scheme, by 
using the input parameters as given in Eqs.(H) and setting \i = 2.5 GeV, we find that: 



d = 1.1245, C 2 
C 5 = 0.0113, C 6 



a 



9 



-0.0096, C 



10 



-0.2662, C 3 = 
-0.0587, C 7 = 
= 0.0026, Cf f 



0.0186, C 4 = -0.0458, 
0.0006, C 8 = 0.0007, 
= 0.3364 



(4) 



where Cf = C 8G + C 5 . 

For the B — > PP decay modes considered here, the decay amplitudes as given in Ref. 
|| will be used without further discussion about details. We focus on estimating the new 
physics effects on those seven well measured decay modes. In the NDR scheme and for 
SU(3)c, the effective Wilson coefficients can be written as || 



a 



i + ^ K + 7v iog- 



a 



a. 



Cj + 2^ A ' i (Ct + C p + C a ) + -£I%C, 



J >j 



71 



(5) 



where A' t = (0, 0, -1, 3, -1, 3, 0, 0, 0, 0) T , B[ = (0, 0, 0, 0, 0, 0, 1, 0, 1, 0) T , the matrices f v and 
jv contain the process- independent contributions from the vertex diagrams [^Q- The 



function C t , C p , and C g describe the contributions arising from the penguin diagrams of the 
current- current Qi,2, the QCD operators Q3-Q6, and the tree-level diagram of the magnetic 
dipole operator Q$g, respectively. The explicit expressions of the functions C t , C p , and C g 
can be found for example in Ref. H . 

In the generalized factorization approach, the effective Wilson coefficients will ap- 
pear in the decay amplitudes in the combinations, 



a>2i 



-1 = u 2?:-i + 



r e ff 
09 



'2i 



n — r e ^ _l_ 



C eff 



2i-l 



5) 



(6) 



where the effective number of colors N^" is treated as a free parameter varying in the range 
of 2 < N%** < 00, in order to model the non-factorizable contribution to the hadronic matrix 
elements. 

In the B rest frame, the branching ratios B(B — > PP) can be written as 



B{B -> XY) = r B 



\P\ 



8vrM| 



\M(B -> XY)f 



(7) 



where t b = 1.653 ps and 1.548 ps for B = B~ and B^ pi[ , respectively, ps is the four- 
momentum of the B meson, Mb = 5.279 GeV is the mass of B u or B^ meson. In the 
numerical calculations we use the same input parameters (the masses, gauge couplings, 





decay constants, form factors, etc.) as in Ref. 
matrix in the Wolfenstein parametrization are: A 



Particularly, the elements of CKM 



0.81, A = 0.2205, 



P 



0.12, 



V 



0.34, 
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which corresponds to 7 = 71° and sin 2(3 = 0.79 favored by the global fit and the new 



measurements 25 



For the seven well measured B — > PP decay modes, currently available measurements 
as reported by CLEO, BaBar and Belle Collaborations [§-f| and their averages are listed in 
Table |. The data have been changed greatly when compared with those in year 2000: 

• For ratio B(B — > 7r + 7r~), the new BaBar result is (4.1 ± 1.0 ± 0.7) x 10~ 6 instead of 
the old (9.3±|:f ±i;l) x 10" 6 . The average therefore decreased to (4.4 ± 0.9) x 10~ 6 . 

• The BaBar measurement of BiB — > i^°7r°) is only ~ 8 x 10~ 6 , the average therefore 
becomes much smaller than two years ago. 

• For the ratio BiB — > K°rj'), both the BaBar and Belle result are much smaller than 
CLEO's measurement, the average is only (56 ± 10) x 10 -6 and clearly smaller than 
the branching ratio of B — > K + i]' decay. 

In the SM and GF approach, the decay B — > -k + ti^ is very simple and receives contribu- 
tions from the dominated tree diagram, the QCD and electroweak penguin diagrams, and 
depend on one form factor F^ n (0) only, as can be seen from the decay amplitude ||: 

M(B° -> 7T+7T-) = -&UFt{ml) [ml - ml 



V2° 

^{v^VZfli-VaVt 



td 



a 4 + a 10 + 2 (a 6 + a 8 ) 



2 

mi 



(m b - m u )(m u + m d ) 



The new physics contribution to this decay in the model III is only 2.5% and thus can 
be neglected. In the SM, there is no other contribution to this decay mode and therefore 
it seems to be a clean decay mode to determine the form factor F^ n (0). For 7 w 71° 
as indicated by the global fit, the interference between the tree and penguin-diagrams is 
constructive, and the measured branching ratio of this decay leads to a small F^ n (0) 

F B7r (0) = 0.21 ± 0.03. (9) 

Which is clearly smaller than the values from Lattice-QCD or Light-cone QCD sum rules: 
F B7r (0) = 0.30 ± 0.04, or 0.28 ± 0.05 as given in Ref. |2§ and Ref. J2?[, respectively. For 



F^ % (0i) = 0.3 and 7 = 71°, however, the SM prediction for the branching ratio is (7.5 — 
10.7) x 10" 6 for N e J f = 2-oo in the GF approach and about 9 x 10~ 6 in the BBNS approach 
|28| , which is clearly too large to be consistent with the data. We thus believe that the form 
factor F S7r (0) should be apparently smaller than 0.3. 

On the other hand, the QCD-penguin-dominated B — > K~n + and K°7r + decays also 
depend on the form factor F s?r (0), as can be seen from the decay amplitudes as given in 
Ref. ||. But one should be very careful to extract the form factor F^^O) from these two 
decay modes because of (a) the neglected rescattering and other non-factorized contributions 
to these two decays may be large; and (b) the new physics contributions in the model III 
are also large, ~ 50% respect to the SM predictions. In the GF approach, the measured 
branching ratios of these two decays prefer a larger F^ n (0): F^lO) = 0.30 ±0.03 if the new 
physics contributions to these two decay modes are not included. But 
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F* n (0) = 0.24 ±0.03. (10) 

if the measured branching ratios of B — > 7r + 7r~, K~tt + , K°tc + decays and the new physics 
enhancements are all taken into account. Although the central value of the form factor 
in Eq.fllOD is still smaller than F B7r (0) = 0.28 ± 0.05 as given in Refs. ||||27|, they are 



compatible within errors. In the numerical calculations, we will use F^ n (0) = 0.24 ± 0.03. 

Furthermore, the form factor F^ K (0) cannot deviate too much from F^ n (0), otherwise 
the SU(3) flavor symmetry will be broken badly. As indicated by the data and theoretical 
considerations, it is a good approximation to take F^ K (0) / F^ n (0) = f k / f T as a measure of 
SU(3) symmetry breaking. We then find that 

F<f*(0) = 0.29 ± 0.04. (11) 

for F Bw (0) = 0.24 ±0.03. 

In the GF approach, the QCD-penguin-dominated B — > Kvl decays depend on the form 
factors Ff K (0) and F W (0), 



F Bv \0) = Ft(0) 



sin 6*8 cos^q 



0.10 ±0.02, (12) 



for F^(0) = 0.24 ± 0.03 , 6 = -9.1° and 6 8 = -22.2° in the two-angle mixing scheme |9 . 



In Table [TT], we show the theoretical predictions for the branching ratios of seven studied 
decay modes. We use the form factors as given in Eqs.( |il3| , |n||T2l) , and keep all other input 



parameters the same as those used in Ref. [|TJj]. The branching ratios collected in Table 
1 1] are the averages of the corresponding B and anti--B meson decay rates. The ratio 5B 
describes the new physics correction on the decay rates and is defined as 

6B(B -+ XY) = b{b - xy)SM (13) 

From Table | and Table y, we find that 

• By using F B7r (0) = 0.24 ±0.03, our predictions for B(B -> tt+tt") in both the SM and 
model III are in agreement with the data. 

• For B — > K7r decays, the SM predictions seem smaller than the measurements. The 
results under the BBNS approach are similar |[22|| . For the B — > K°n decay, specifi- 



cally, the SM prediction is about half of the measured decay rate, and the new physics 
enhancement is essential for the theoretical prediction to become consistent with the 
data. 

By taking into account the uncertainties of those input parameters as given explicitly 



in Ref. [14], we find numerically that 



where the uncertainties of those input parameters have been considered. It is easy to 
see that the SM predictions in the GF approach is about half of the measured value. 
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In BBNS approach, the theoretical predictions for the corresponding branching ratios 
are also much smaller than the experimental data . The new physics enhancement 
can boost the theoretical predictions close to the lower part of the measured values, 
but still leaves a moderate space for additional contributions. 

Since 1997, the unexpected large rf production has been widely discussed in the literature 
PH|- P3|PU| - P^|J2E[ | . For the sake of completeness, we here make a brief comment on some 



typical interpretations in the framework of the SM. 



Atwood and Soni |IJ gave arguments for the need of enhanced b — > sg* decays followed 



by g* — > rfg via the QCD gluon anomaly. Taking a constant ggr/' vertex form factor 
i7(0, 0, 77?y), the observed large branching ratio B — > 1]'X S can be achieved. But as 
pointed by Hou and Tseng fH|, Kagan and Petrov ||15|| , if one considers the running 
of a s , and the m?,/(q 2 — ml,) dependence of the ggrf coupling, the result presented in 
PHI will be reduced greatly. 

Halperin and Zhitnitsky [T3| argued that the dominant contribution to r/ production 
is due to the Cabibbo favored b — > (cc)is process followed by the transition (cc)i — > rf, 
i.e. the "intrinsic charm" component of rf . But according to the explicit calculations in 
Refs. this mechanism can not give a good explanation for the measured B — > Krj' 
decay rates. 



Yuan and Chao [IT| argued that the inclusive rf production in B decays may dom- 
inantly come from the Cabibbo favored b — > (cc)gs process where cc pair is in a 
color-octet configuration, and followed by the nonperturbative transition (cc)§ — > rj'X. 



• The authors of Ref. (Tj| proposed di-gluon fusion mechanism. It seems that this 
mechanism could enhance 7/ production. But because of our ignorance about the form 
factor of g*g*T]' vertex, there are large uncertainties in calculation. 

• Recently, M.Z.Yang and Y.D.Yang reconsidered the di-gluon mechanism and gave a 
calculation for B — > rf'P in the BBNS approach. They computed the vertex g*g*rj' 
in perturbative QCD, and found that the branching ratios of B — > Krj' are really 
enhanced and in agreement with the data. But, as pointed in Ref. [p^] , the consistency 
of Yang's perturbative calculation is questionable due to the end point behavior. 

From above discussions, one can understand that it is still an open question for us to interpret 
the large branching ratios of B —>■ Krf decays. Further investigations for various possible 
mechanisms are welcome. 

Furthermore, because of the isospin symmetry between the u and d quarks, the decays 
B — > K + f]' and K°r] f should theoretically have similar branching ratios. The known new 
mechanisms in the SM or the new physics models also contribute to these two decay modes 
in very similar way. The measurements of the BaBar and Belle collaboration, however, 
show a clear difference between these two decay rates. We do not know how to interpret 
this difference. One may need something new to resolve this problem if it is confirmed by 
the forthcoming data. 

In short, we reexamined the branching fractions of seven well measured B — ► PP decay 
modes in the SM and model III, and compared the theoretical predictions with the new 
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data. Within the considered parameter space we found that: (a) the measurements of the 
branching ratios for B — > 7r + 7r~, K~tt + and K°tt + lead to a strong constraint on the form 
factor F^fo): F^ n (0) = 0.24 ± 0.03; (b) the new physics enhancements to the penguin- 
dominated B — > Kir and Krj' decays are significant in size, ~ (40 — 65)% w.r.t the SM 
predictions, and play an important role in restoring the consistency between the data and 
the theory. 
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TABLES 



TABLE I. Measurements of the branching ratio B(B — ► PP) (in units of 10 6 ) as reported by 
CLEO, BaBar and Belle Collaborations. The last column lists the world average. 



Decay Mode 


CLEO 


BaBar 


Belle 


Average 


B° tt+it- 


4.3l};jj ± 0.5 


4.1 ±1.0 ±0.7 


5.6t£j| ± 0.4 


4.4 ±0.9 


B° -► ivT+vr- 
B + tf+7r° 

#o ^ ^0^0 


17.2^24 ± 1-2 
11 6 +3:o+i.4 

ll.D_2 7 _! 3 

18.2+^ ± 1.6 


16.7 ± 1.6 ± 1.3 
10.81^9 ± L0 
18.2l|;g ± 2.0 
8.2l| 7 ± 1.2 


1Q 0+3.4+1.5 

iy - o -3.2-0.6 

i6.3iflirJ 

10 7 +5.7+1.9 
J-" 3 -' -4.8-1.8 
1fi n+7.2+2.5 
iD - u -5.9-2.7 


17.3 ± 1.5 
12.1 ± 1.7 

17.4 ±2.6 
10.7 ±2.7 


5+ if+r/ 


80i£° ± 7 
89+{^ ± 9 


70 ± 8 ± 5 
42+^ ± 4 


79+^ ± 9 

551" ± 8 


75 ±7 
56 ± 10 



TABLE II. Branching ratios (in units of 10 6 ) of seven B — > PP decay modes in the SM and 
Model III by using F B7r (0) = 0.21 (the first entries), 0.24 (the second entries), 0.27 (the third 
entries), and assuming iV c e // = 2 - 00, 6 = 0° and M H + = 200 GeV. The last column lists the 
world average data. 







SM 






Model III 




SB [%} 


Data 


Decay Mode 


2 


3 


00 


2 


3 


00 


N^t = 3 




B° -► vr+vr- 


3.66 


4.16 


5.24 


3.75 


4.26 


5.37 


2.5 


4.4 ±0.9 




4.78 


5.43 


6.85 


4.90 


5.56 


7.01 


2.5 






6.05 


6.87 


8.66 


6.20 


7.04 


8.87 


2.5 




B + _> K+tt 


4.93 


5.53 


6.85 


7.12 


8.01 


9.98 


44.9 


12.1 ± 1.7 




6.44 


7.22 


8.95 


9.29 


10.5 


13.0 


44.9 






8.15 


9.14 


11.3 


11.8 


13.2 


16.5 


44.9 






7.18 


7.96 


9.66 


10.8 


12.0 


14.7 


51.3 


17.3 ± 1.3 




9.37 


10.4 


12.6 


14.1 


15.7 


19.2 


51.3 






11.9 


13.2 


16.0 


17.9 


19.9 


24.3 


51.3 




5 + 


8.12 


9.46 


12.5 


12.1 


14.1 


18.4 


48.9 


17.4 ± 2.6 




10.6 


12.4 


16.3 


15.9 


18.4 


24.1 


48.9 






13.4 


15.6 


20.6 


20.1 


23.3 


30.5 


48.9 




5° - K 7T° 


2.90 


3.32 


4.25 


4.56 


5.21 


6.66 


56.9 


10.7 ±2.7 




3.80 


4.34 


5.56 


5.96 


6.80 


8.69 


56.9 






4.81 


5.49 


7.03 


7.54 


8.61 


11.0 


56.9 




B+ X+r/' 


9.69 


12.2 


18.2 


16.3 


20.1 


29.0 


64.7 


75 ±7 




12.7 


16.0 


23.8 


21.3 


26.3 


37.9 


64.7 






16.0 


20.2 


30.5 


27.0 


33.3 


47.9 


64.7 






9.33 


12.0 


18.3 


15.6 


19.5 


30.0 


62.8 


56 ± 10 




12.2 


15.6 


23.9 


20.4 


25.5 


37.4 


62.8 






15.4 


19.8 


30.2 


25.8 


32.2 


47.3 


62.8 
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